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Abstract 

We present a simple construction that maps quantum circuits to graphs and vice-versa. 
Inspired by the results of D.A. Lidar linking the Ising partition function with quadratically 
signed weight enumerators (QWGTs), we also present a BQP-complete problem for the 
additive approximation of a function over hypergraphs related to the generating function 
of Eulerian subgraphs for ordinary graphs. We discuss connections with the Ising partition 
function. 
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I. INTRODUCTION 



Relationships between quantum computation and graph theory are emerging beyond ap- 
phcations to graph theoretic problems. For example in the one-way quantum computation 
setting, quantum graph states and their relationship to entanglement is already well known 
|15j . Here one sees the correspondence between a graph and a quantum state where the 
vertices correspond to the qubits of the state and the edges correspond to pairs of qubits. 
For quantum circuits one may make use of a graph to actually represent the circuit or ar- 
chitecture [9]. For example, if we let F = {V, E) be a graph, then the set of vertices V may 
represent the individual qubits (or input into the circuit) and the edges E correspond to any 
pair of qubits that may be acted upon by a two qubit gate. In another related approach 
provided in [23], they instruct to "regard each gate as a vertex, and for each input/output 
wire add a new vertex to the open edge of the wire." Using this correspondence they prove 
statements about families of quantum circuits that are classically simulatable. For the work 
here we assume the usual quantum circuit model and that one may efficiently swap between 
different universal gate sets due to the famous theorem of Solovay and Kitaev [12]. The 
methods used here have been employed to study quantum circuits which may be simulated 
classically [1]. 

In this paper we provide a simple construction which makes a direct connection between 
quantum circuits and graphs via their incidence structure. As an application we construct a 
function related to the generating function of Eulerian subgraphs and discuss the relationship 
to the Ising partition function. We show that quantum computers can provide additive 
approximations of this related generating function, which we call the signed generating 
function of Eulerian subgraphs, E'{r, A). We demonstrate that it is a BQP-complete problem 
(when we allow it to be defined over hypergraphs which are a generalization of graphs) as 
it is intimately related to quadratically signed weight enumerators (QWGTs) [16j via this 
construction. It is well known that the Ising partition function Z may be expressed in terms 
of the generating function of Eulerian subgraphs, E{G,X). We provide some ideas for the 
future to use -E"(r, A) for efficient additive approximations of Z. Recently in [7|, an additive 
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approximation algorithm of the Tutte polynomial was given which solved instances shown 
to be BQP-complete. As the Ising partition function is just a speciahzation of the Tutte 

polynomial, the complexity of additive approximations of certain non-planar instances of the 
Ising partition function is an interesting open problem. 

A. Ising spin model 

We shall first introduce the Ising spin model in order to motivate the computation of gen- 
erating function of Eulerian subgraphs. Let G = {E, V) be a finite, arbitrary undirected 
graph with \E\ edges and \V\ vertices. In the Ising model, each vertex is associated with a 
classical spin (ctj = ±1) and each edge G E with a bond {Jij — ±J). The Hamiltonian 
of the spin system is 

The probability of the spin configuration a in thermal equilibrium at temperature T is 
given by the Gibbs distribution: P{a) = ^W{a), where the Boltzmann weight is W{a) = 
exp[—(3H{a)], (3 = l/kT, and Z is the partition function: 

cr 

B. Generating function of Eulerian subgraphs 

An Eulerian subgraph of a graph F is a set of edges that form a tour (a path that begins and 
ends at the same vertex) in which every vertex is of even degree. The generating function of 
Eulerian subgraphs of T is given by 

E{T,x) = 

a 

where the sum is over all Eulerian subgraphs and wt is a weight function (in this case the 
number of edges in the subgraph) . This brings us to another expression for the Ising partition 
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function which was discovered by van der Waerden and is given by 

Z(/3) = 2l^l Y[ cosh{/3Jij)E{T,tanh{(3Jij)) (3) 

It is this form of the partition function that motivates this work [221 126] . 

II. QWGTS AND THEIR RELATION TO THE ISING PARTITION FUNCTION 

Definition 1 Quadratically Signed Weight Enumerators (QWGTs) in general are of the 
form fWj 

SiA,B,x,y)= J2 i-ir'x^S''-^^ (4) 

b:Ab=0 

where A and B are 0, 1-matrices with B of dimension nxn and A of dimension mxn. The 
variable b in the summand ranges over 0, 1-column vectors of dimension n. All calculations 
involving A,B or b are done modulo 2. 

Note that the evaluation of a QWGT, given x and y are natural numbers, is in general a 
T^P problem and as it includes evaluations of the weight enumerator of binary linear codes 
it is in fact #P-complete [12]. 

We shall now review how QWGT's were constructed in [TB] in some detail. Let G be a 
quantum circuit and U (G) the corresponding unitary operator. Note that a universal gate 
set can be achieved by allowing arbitrary rotations about any product of Pauli operators i.e. 

e-^-^/^ = cos (^)^ + ^sin(^0 a, 

where ab = nr=i "^bf ^'^'^^ that ctoo = I,<^oi = cm = cry and aio = oz [IE]. This means 
that 6 is a binary vector whose length is 2n, twice the number of qubits, and the superscript 
ii) represents the qubit which is operated on by the corresponding Pauli matrix. 

It is possible to express our unitary operator as a product of real gates and we can do this 
as follows. Take the product U (G) = GnGn-i ■ ■ - Gi where each gate is of the form 

1 

Gfc = - (a ± i/^o-foj 
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where again bk is a binary vector of length 2n but each bk can only contain an odd number 
of ll's i.e. each gate can only have an odd number of Pauli Y operators ay- Also note that 
7 = Y^a^ + As a further modification, which will allow us to have simple multiplication 
rules for our gates, define 

where |6|y is the number of cry's occurring in a^^. 
We can now write 

7 

Now define C to be the block diagonal matrix whose blocks consist of 




Then the property that has an odd number of ll's is given by b^Cb — 1 and so we have 
the multiplication rule 



= {-lf^^''ak,+b, (5) 
where the addition in the subscript is bit by bit modulo 2. 

Let H be the {2nxN) matrix whose columns are the bk- H isa. polynomial size representation 
of the quantum circuit where each column represents a gate and every pair of rows represents 
a qubit. We then have the following expansion. 

1 

U{G) = n (6) 

k=N 

1 1 

k=N ' 

i_ ^(_l)a*lwtr(i/*C//)a^|a|^iV-|a| (g) 
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Now note that if we only sum over the a's such that CHa = then we assure that 

(00 ■ ■ ■ |f/(G)|00 ■ • ■ ) = ^ ^(_i)a'lwtr(//*CH)a^|a|^Ar-|a| 

^ a 

is always non-zero for this omits X and Y gates from our sum. 

As a simple example to illustrate the correspondence between the matrix representation H 
of the circuit and the actual operation of the circuit, consider 



H 



1 1 1 

1 

1 1 

1 

1 1 1 

1 1 

This matrix represents a circuit which operates in the following way: 

g-iz(i)(g.x(2)(g.y{3) I g-jz(i)(g.z{2)(g.y(3) | g-iy(i)(g.z(2)(g.z(3) f 

where the superscripts represent which qubit is being acted upon. Thus each column encodes 
each exponentiated operator, i.e., each gate. When using our proposed gate set, we would 
have 



^[{al - IpZ^^^ (g) X(2) ® _ ^^^(1) ^ ^(2) ^ y{3))(^j _ -pyW ^ ^(2) ^ ^(3))] 

In [T^ it was shown that the Ising partition function can be expressed in terms of a QWGT. 
Let A be the incidence matrix of a graph g, i.e., 

. 1 (v = i and (i, ?') G E) 
else 



Then we have 

o\v\ o\v\ 

^ ^ a ^ ^ aekcxA 

S{A,dg{w),X,l) (11) 



;i - A2)i^i/2' 
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where w = (wi2, W13, . . . ) {w gives the distribution of ferromagnetic {wij = 0) or anti- 
ferromagnetic {wij = 1) interactions along the edges of the given graph), A = tanh(/3J) (the 
"temperature"), V is the set of vertices, E is the set of edges and B = dg{w) is the diagonal 
matrix formed by putting w on the diagonal and zeros everywhere else. This form of Z will 
be considered in future work. 



III. A RELATIONSHIP BETWEEN HYPERGRAPHS AND QUANTUM 

CIRCUITS VIA QWGTS 



The following "mapping" between hypergraphs and quantum circuits was first introduced in 
order to find a way to compute the Ising partition function. It was extended to obtain a class 



of quantum circuits which can be simulated classically in |1] and is based on equation (11). 
This mapping involves interpreting the matrix representation H of the quantum circuit, as 
outlined above, as coming from the incidence matrix of a hypergraph. In this way we have 
a many to one mapping from quantum circuits to a hypergraph. 

First we define hypergraphs. 



Definition 2 A hypergraph is a generalization of a graph where edges are replaced by hy- 
peredges. Let V = {vi, t>2, . . . , Vk} be the set of vertices and let E = {61,62, ■ ■ ■ , 6n} be the 
set of hyperedges. Each 6i = {va, Vi2, ■ ■ ■ , Vim} is a collection of vertices where each Vij G V . 

The standard reference for hypergraphs is [30J. Note that graphs are just a special case of 
hypergraphs where each edge just consists of two vertices and that the incidence matrix is 
defined in the same manner as above. 

For the description of the mapping in the next subsection, we shall restrict our quantum 
circuits so that the corresponding hypergraphs are ordinary graphs. We shall point out when 
this restriction is important. 
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A. THE MAPPING 



The motivation for this mapping is to obtain a QWGT equal to a matrix element of the 
unitary matrix of a quantum circuit that looks something like the generating function of 
Eulerian subgraphs E{r, A). If we were able to efficiently approximate -^(r, A) then according 
to equation ^ we would have an efficient method of approximating the Ising partition 
function. 

It turns out that if we take the ansatz 



G. = ^^=^(A + (12) 



for the gate set we obtain 



1 1 



k=N 



^Il(A + ^6j (14) 



(A^ + 1) 



k=N 



(A2 + 1)^/ 

Now, ignoring the normalization we have 



1 ^(_,).nwM^^CH).^H~^^_ (15) 



(00---|;7(G)|00---) oc (-l)'^*i"*''(^*^^)'^Al'^l 



a(^Ker{CH) 

Let us make the following assumptions: 



1. Take CH to be a binary matrix with only two I's per column i.e. CH will be identified 
with the incidence matrix of the given graph F. For hypergraphs, this is not necessary 
as the columns of a hypergraph incidence matrix may be populated by more than two 
1 's, as this corresponds to edges consisting of multiple vertices. 

2. is a matrix of dimension 2n x N with one (11) and at most one (01) per column, 
i.e., one Y operation and at most one X operation per gate respectively. The source 
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of this restriction will be explained below. For example a column may look like 

( 110001100010 ) . 

H encodes the quantum circuit. For hypergraphs, these restriction vanish, however it 
is necessary that there are an odd number of (11) 's per column as our gate set depends 
on this restriction. 

These assumptions will provide the basis for a natural mapping between quantum 
circuits and graphs. Consider one additional assumption. 

3. 

a*lwtr(//*C//)a = mod 2 Va e ker(C//). (16) 

This ensures that the matrix element (00 • • • |?7(G) |00 ■ ■ ■ ) is equal to the generating function 
of eulerian subgraphs. This is achieved as follows. 

First we need to associate the incidence matrix A with a matrix CH. The only thing we need 
to do is to create a matrix with double the number of rows of A, with row 2i — 1 occupied 
by the i*^ row of A, and each even row the zero vector. Thus we obtain 



/ 


An 


A12 .. 


• AiN 


\ 







.. 


. 






A21 


A22 ■■ 


■ A2N 




V 





.. 


. 


/ 



As far as the graph is concerned, this amounts to adding isolated vertices, which does not 
add any cycles. We now have the 2n x N matrix CH as our representation for P. By the 
action of C, we see that CH gives us some freedom in our choice of H, which is the matrix 
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representation of the quantum circuit. Specifically we have 



Xi 


X2 . 


■ Xn 


An 


A12 . 


■ AiN 


Xn+1 


Xn+2 ■ 


■ X2N 


A21 


A22 . 


■ A2N 



\ 



H 



\ Ani An2 ■ ■ ■ Ann ) 

The Xi must be selected according to the constraints mentioned above. We see that column 
k oi H will only have two Aii^'s which are equal to 1 as these come from an incidence matrix 
and by definition, an incidence matrix has only two I's per column. One 1 is possible as this 
represents a loop, i.e., an edge that begins and terminates at the same vertex. Further, by 
the QWGT formalism constructed above, we must have an odd number of Y ((11) entry in 
the column) operations per gate. Hence, there must be one 11 per column in the matrix H. 
There are only two positions where we could select an Xi in column /c to be 1 such that it 
will be followed by an Ajk that is equal to 1. Thus each column (or gate) must have only 
one Y operation. By the same reasoning we see that there is only one possible place to put 
an X operation, i.e., only one way to place a (01) in column k. So there can be at most one 
X operation per gate. There is no restriction as to the number of Z operations per gate as 
we have the freedom of putting a 1 before any Aik that is set to 0. By turning the Xi on or 
off (1 or respectively) we obtain different circuits. This provides a degree of freedom that 
allows one to choose a quantum circuit that may satisfy the final assumption which ensures 

that the sum EaeKer(CiJ)(-l)"*'"*'^'^*^'^^"'^'"' is equal to EaeKer(Cif) which is E{r,X) as 
desired. 

Without the restriction given by 

a*lwtr(i/*Ci/)a = mod 2 Va G ker(Ci/), 

we actually have that 

(00---|t/(G)|00---) 



(A2 + 1)1^1/2 



J2 (-l)''"Al"l. 



aeKei{CH) 
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This means that knowledge of the matrix element (00 ■ • • |f/(G)|00 ■ ■ ■ ) amounts to knowledge 



where the ha refer to the a*hNii{H^CH)a. We shall call E'{T ^ A) the signed generating 
function of Eulerian subgraphs as the sum is over all subgraphs whose edges have even 
degree. Specifically, as CH is associated with the incidence matrix of a graph, the whole 
null space of CH are the characteristic vectors of all Eulerian subgraphs |28j. 



Theorem 1 Additive approximations of E'{r,X) over hypergraphs is BQP-complete. 

Before proving this theorem we would like to note that we include evaluations of E'{r,X) 
for hypergraphs F only for completion, but that this is not necessary. The restriction to 
graphs forces the corresponding gates to allow an X operation on one qubit, and forces one 
to have a Y operation on another qubit , but an arbitrary number of Z operations on the 
remaining qubits. As these gates correspond to exponentiated Pauli operators, these are 
multi-qubit operations and thus it is easy to implement entanglement under this restriction 
as well as control gates. Thus, from the results in [311 [32] we see that the quantum circuits 
corresponding to ordinary graphs are capable of universal quantum computation. In addition, 
as our mapping depends on the sum over all simple cycles of a given graph, any one qubit 
operation may be inserted without effecting the sum, as these correspond to adding loops, 
i.e., an edge that begins and ends at the same vertex. 

Proof of Theorem 1: 

It is clear from the above that an approximation of the matrix element (00 ■ ■ ■ \U (G)|00 • ■ ■ ) 
will give an approximation to E'(r,X). Recall from [HI [2^ that via the Hadamard test 
one can obtain an additive approximation of this matrix element. This means that using a 
quantum computer we can retrieve an approximation m (with some probability of success 



of 




(17) 



aeKeriCH) 



IV. BQP-COMPLETENESS 
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bounded below by .75, say) which satisfies 

(00 ■ ■ ■ |f/(G) |00 ■ ■ ■ ) - p < m < (00 • • ■ |t/(G) |00 ■■■)+p 

where p is a polynomially small parameter. (This intuitive form can be easily derived via 
the definition of an additive approximation provided in [1].) 

Now, assuming that we have at our disposal the universal gate set given by ^ = 
±2 arccos(4/5) rotations of products of Pauli operators, then with an overhead of 
polylog(|-E|/e) of gates we may approximate our gates Gk, to accuracy 0{e/\E\) [TB]. This 
means that we may indeed approximate the signed generator function of Eulerian subgraphs 
via the Hadamard test and so this problem is in BQP. 

Now, the other direction. We must demonstrate that knowledge of -E"(r, A) is enough to 
simulate any quantum circuit. First, any quantum circuit corresponds to a hypergraph 
under the scheme presented above. Since BQP is a decision class all we have to do is convert 
a quantum circuit into its decision making counterpart. This may be done as follows |21j . 




10, 



FIG. 1: A circuit illustrating the procedure to apply a quantum circuit [/ to a decision problem. 

Let U he a quantum circuit for some decision problem and without loss of generality, assume 
that the "yes" or "no" answer is given by the output of the first qubit Ui, i.e., |0) or |1) is 



13 



"yes" or "no" respectively. The remaining qubits are ignored and assumed to be extraneous. 
Now, take an ancilla qubit Ua set to |0) and adjoin it to U. Next, CNOT the output of Ui 
with Ua and encode the answer as the state Apply the inverse of the circuit, U'' , to all 
the output qubits except lip), to uncompute the outputs of U to |00 ■ ■ ■ 0). In this way, one 
arrives at the state \00---0)\4') which will either be |00---0)|0) or |00---0)|1), effectively 
"deciding" the decision problem. Thus one can assume, with no loss of generality, that any 
quantum circuit that solves some decision problem either outputs |00 • • • 0)|0) or |00 ■ • • 0)|1) 
(see FIG.[T]). This argument depends on the quantum circuit being able to output the correct 
answer with certainty. This is of no matter as a similar argument can be made for a circuit 
which outputs the correct answer with some constant probability above a half [27] . Thus, 
knowledge that (00 ■ ■ ■ |[/|00 ■ ■ ■ ) = 1 implies that the control qubit will be |0) and thus 

1^) = |o). 

This means that knowledge of -E"(r, A) can be used to effectively decide the decision problem 
for any quantum circuit as it is proportional to the matrix element (00 ■ ■ ■ |?7|00 ■ ■ ■ ). In our 
case, a natural decision problem would be to decide if E'(T, A) is bounded above by some 
constant or to decide its sign. The result of this decision would correspond to {tp) either 
being |0) or |1) in FIG.[i|B 

Important Caveat: Upon careful inspection of the above argument one may find something 
amiss. We are referring to the idea that perhaps knowledge of E'{r,X) may in fact not be 
enough to solve all decision problems in BQP, since only one A is specified. Recall that our 
universal gate set consists of rotations about products of Pauli operators, with two angles 



at our disposal, namely ±2 arccos(4/5). When we use the gate set given by the ansatz (12) 



the temperature A plays the role of the angle, and so it seems that there is only one angle 



available for the gate set Gk- In fact, gates of the form (12), consist of multiple angles [16] . 
A quick calculation verifies the following claims. 

1. If one choses A = |, then one indeed recovers the rotational angles 9 = ±2 arccos(4/5). 
However, in the scheme outlined in this paper this is not acceptable when one moves 
to the Ising model, as this particular choice of A means that the physical quantity 
l3J become complex. In fact one would need j3J = l2i^±H!:_ This is of no matter for 
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the quantity E'{r,X) but is not acceptable when considering its interpretation as a 
partition function. However, this does mean that the gate set given by our ansatz is 
capable of universal quantum computation |16j . 

2. If one chooses, for example, A = | then one has access to the two rotational angles 
6 = ±2 arcsin(4/5) and thus, we may indeed claim universality for our gate set given 



by equation (12). Further, all quantities are now physically acceptable. 



A. Examples 



Here are two simple yet instructive examples. 
1) Let the given quantum circuit be encoded by 



H = 



1 

















1 








1 











1 

















1 








1 
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1 


1 








1 
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1 











1 
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1 


1 


1 









The incidence matrix (ignoring isolated vertices) can be retrieved easily from H and it is 
given by 



^ 1 1 0^ 

10 10 

10 11 

^ 1 1 1 1 Oy 

The corresponding graph is given in FIG. [2j 
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FIG. 2: The graph obtained from the simple circuit given by H. 

2) A simple demonstration of a controlled 2-qubit gate is given by the sign flip operator [1 
given by 



which acts in the following way 



^2 7r , 



1 06) — > e-*"-?|06) 
and 

1 16) —^e'''"^lh). 

In our gate set we would have the corresponding gate, al — i[3a\ ® a1 with the action 

|06) — > (a-i/?)|06) 

and 

|16) (a + z/3)|16). 
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V. FUTURE WORK: APPROXIMATING THE ISING PARTITION FUNCTION 

Recall the following form of the Ising partition function 

Z(/3) = 2l^l Yl cosh{l3Jij)E{T,t&nh{f3Jij)) 

and note the difference between the function we are able to approximate via quantum com- 
putation, E'{T,X), and the actual generating function of Eulerian subgraphs E(T,X). We 
have that 



E'{T,\)= ^ l^_-^yHwtv{H'CH)a^\a\^ 
aeKer(CH) 

If we wanted to use this for the approximation of the Ising partition function, as previously 
mentioned we would require that a*lwtr(/7*Ci/)a = mod 2 Va G kei^CH). If this 
requirement was met then we would run the quantum approximation algorithm with gate 
sets corresponding to A = tanh(/?Jjj)). This would require Od-El) different approximations 
as indicated by the product over all edges in the above formula. In effect we would have 
a polynomial additive approximation of the Ising partition function for any set of edge 
interactions and for any graph. 



But alas, there is a problem. The equation that must be solved (equation (16)) in order 
to ensure that E'{r, A) = E{r, A) in fact determines which particular quantum circuit must 
be used for the computation. By brute force this could require an exponential number of 
calculations in the number of vertices. Future work will involve studying this approach to 



see if one can in fact a priori guarantee that equation ( 16 ) is satisfied for certain non-planar 
graphs. For example, if one has knowledge about the parity of all the Eulerian subgraphs 
(number of edges) then this may be used to efficiently find the representation H of the 
quantum circuit required. The cubic lattice is an example where every Eulerian subgraph 
contains an even number of edges. This issue also arises in a very similar approach outlined 



in |1] but which deals explicitly with equation (11). Other applications of E'{r,X) will be 
explored as well as an extension to a two variable function. We will also attempt to use the 
methods here to find the instances of the Ising model for which evaluations of the partition 
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function is BQP-complete. 



VI. CONCLUSION 



We provide a new way of relating quantum circuits to graphs and vice- versa via an incidence 
structure of the circuit or graph. We also provide a generating function related to the 
generating function of Eulerian subgraphs and demonstrate that additive approximations of 
it for hypergraphs are BQP-complete. Connections to the Ising spin glass partition function 
were made and a discussion of future work dealing with additive approximations of the Ising 
partition function was provided. 
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